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Abstract 

The process e + e~ — > W + W _ — > fif 2 4 contains charges in the initial, 
intermediate and final stages. This gives a rich selection of possible QED 
interconnection effects. Coulomb interaction is the simplest of these, and 
can thus be used to explore consequences. We study a number of experi- 
mental observables, with emphasis on those related to the W momentum 
distribution. Second-order Coulomb effects are shown to be practically neg- 
ligible. The limited LEP 2 statistics will not allow detailed tests, so any 
theory uncertainty will be reflected in the systematic error on the W mass. 
Currently the uncertainty from this source may be as high as 20 MeV. 



1 Introduction 



A precision measurement of the W-boson mass is one of the main objectives of the 
LEP 2 e + e~ collider. To exceed the precision of pp collider experiments, Mw should 
be measured to an accuracy of 50 MeV or better, see e.g. ref. The success of such 
an ambitious goal relies on an accurate theoretical knowledge of the dynamics of the 
production and decay stages in e + e~ — > W + W~ — > 4 fermions. Moreover, owing to 
the large W width T w , these stages are not independent but may be interconnected by 
QED and QCD interference effects, which also must be kept under precise theoretical 
control. Interconnection phenomena may efface the separate identities of the two W 
bosons, and the final state may no longer be considered as a superposition of two separate 
W decays. All such effects related to the radiatively corrected (total and differential) 
W + W~ production cross section should be accurately calculated, so that the W mass can 
be extracted from data. 

For the totally inclusive W + W~ production cross section there is a general proof 
|2|. |3| (also H) that the radiative interconnection effects are suppressed by 0(aT w /M w ) 
or 0(a^T^/M^). The only exception is the Coulomb interaction between two slowly 
moving W's, which is modified by the instability effects in the very narrow region just 
near threshold, \fs — 2M W ^ T w . By contrast, differential distributions could be distorted 
on the level of 0(a) or 0(a1). (Analogous QCD final-state interaction effects appear in 
the differential distributions of tt production processes |3], £|, |, ||, [7].) These distortions 
should vanish in the zero- width limit. In ref. || we studied the possible impact of 
QCD interconnection (colour rearrangement) effects for the direct reconstruction of the 
W mass from the hadronic decay products (as favoured by experimenters |§). The limited 
understanding on the hadronization of overlapping hadronic systems was shown to imply 
a systematic error on of 30-40 MeV. Additional errors could come from Bose-Einstein 
effects among the final-state hadrons [pi], [Tl| . 



There are also purely QED final-state interactions, however, that can produce non- 
factorizable radiative corrections to the Born cross section of e + e~ — > W + W~ — > 



4 fermions [12, El. Specifically, there is a class of virtual interference effects correspond- 



ing to the so-called charged-particle poles 0, 0L O , which induce an explicit dependence 



of the differential cross section on the W-boson virtualities. This is related to the fact 
that the dominant contribution to the non-factorizable radiative interferences comes from 
sufficiently soft photons with k 1 ^ T^. The phenomenon could be exemplified by the 
Coulomb interaction between two unstable W bosons |2|, [l4j] (also fll5l). (To the best 



of our knowledge, the necessity to take into account the effects of Coulomb final-state 
interactions on the W momentum distribution was first pointed out in ref. |§.) Coulomb 
terms could be uniquely separated from the other electroweak corrections in the threshold 
region, where they provide the dominant source of the off-shellness effects. 

Other radiative mechanisms may become essential at larger W-boson velocities /3, 
for example, those caused by the intermediate-final or final-final interferences |3], [13| . 



(We recall, however, that the Coulomb correction is not uniquely defined and gauge 
independent in the relativistic region.) In the ultrarelativistic region it is expected 
that the dependence on the W virtuality disappears in the full expression for radiative 
corrections. This phenomenon is deeply rooted in the conservation of a charged current. 

In the threshold region the Born cross section depends strongly on the momentum p 
of the unstable particle. This leads to a competition between the Breit-Wigner factors 
and the rest of the cross section (oc (3) M, which induces a significant change between the 
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nominal and the actual average mass. In the energy range 

m w > e > r w , (i) 

where E = yfs — 2M\y is the nonrelativistic energy of the two W's, the mass shift towards 
lower values is of 0{T%j/E). This shift changes sign at around 190 GeV because of 
the increasing role of the negative higher-order terms (oc (3 2 ) in the Born cross-section 
expression [[RJ. 



Our interest is in the dynamical QED effects, on top of the purely kinematical dis- 
tortions. The Coulomb final-state interaction between unstable W's could, in principle, 
induce a systematic shift in the reconstructed Mw of O(craTw) [14], [31| that is of the 
same magnitude as the aimed-for precision of LEP 2 measurements. Therefore the shifts 
in the W momentum distribution is the result of a complicated interplay between phase- 
space and virtuality-induced effects appearing in the radiatively corrected cross section 
for e + e~ — > W + W~ — > 4 fermions. In the main region of LEP 2 running, y/s — 175 GeV, 
the two effects above are of the same order. To obtain reliable numerical results one also 
needs to analyze carefully the role of other final-state interaction mechanisms (for such 
an attempt see ref. |13|). 



In this paper we analyze the impact of QED interconnection, exemplified by the 
Coulomb effects on the W momentum distribution in the threshold region. This region 
is relevant for LEP 2 physics, and it is just here that off-shell and finite- width effects are 



most important [0, BL 14, ujn. It may be seen as a continuation of the studies presented 



for mass distributions in ref. 1T6 



Particular emphasis is put on the momentum distribution of the W's. This is especially 
relevant since the Coulomb effects depend primarily on the W momentum. However, there 
is also an experimental reason for such a choice. In the best of possible worlds, the final 
state could be subdivided into two sets of particles, one stemming from the W + decay and 
the other from the W - . The invariant masses mi and m 2 could then be calculated and 
directly compared with theory. In addition to the assignment problems already mentioned, 
such an approach suffers from the limited energy resolution of detectors. Jet directions, 
on the other hand, are much better reconstructed. If energy and momentum conservation 
relations can be used, i.e. if jet energies are suitably rescaled to the known e + e~ beam 
energy, the W mass is better constrained ||. Optimal use is made of the experimental 
information if the W + and W~ of each event are assumed to have the same mass, i.e. if one 
single average mass m is extracted for each event, rather than m\ and m 2 separately. Such 
a mass measurement obviously is equivalent to a momentum measurement, according to 



the relation y/s = 2\Jm 2 + p 2 . However, the momentum is the more direct observable, 
since it is related to the acollinearity of the jet pairs; we recall that the limit p = 
corresponds to each hadronic W decay giving a pair of back-to-back jets. A measurement 
along these lines is therefore excluded only when both W's decay leptonically. 

The results presented here are aimed at the rather modest task of evaluating the size 
of the QED interconnection effects. They should only be considered as a qualitiative guide 
rather than a complete numerical prediction. In what we follows, incoming electrons and 
positrons are assumed unpolarized, and initial-state radiation (ISR) effects are neglected. 
We also assume that all particles can be perfectly well measured and neglect the QCD 
interconnection effects. 

In section 2 some of the basic Coulomb formulae are reviewed, with emphasis on 
the physical origin and consequences. Section 3 contains a numerical study of several 
observables relevant at LEP 2, with conclusions on the practical importance of Coulomb 
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effects, on the impact of theoretical uncertainties, and on the choice of experimental 
observables. Finally, section 4 contains a summary and an outlook. 



2 Qualitative discussion 

For illustrative purposes we present here the results for momentum distributions in the 
non-relativistic approximation. The numerical calculations in the next section, however, 



are based on the full relativistic formulae of ref. [14 



As shown in ref. |7j (see also ref. (T4J) the non-relativistic Green's function describing 
the interaction of the W's depends only on one off-shell variable, namely the W momentum 
P (p — |p + l = |p~| i n the rest frame of the pair). The integration over the W squared 
masses s\ and S2 reduces to an integration over p 2 with 



V 



[s-s 1 - s 2 ) 2 - 4gjg 2 
4s 
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In the dominant p region the differential distribution for e + e 
be written in terms of the rapidly changing variable 



W+W" 
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M W E - p 2 
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(2) 

fif 2 f 3 f 4 can 
(3) 
(4) 



with p = Ipjyfs. Here 5u is the so-called hard correction and a5c/P is the first-order 
Coulomb term B Ol , 



5c 



7T \k\ —p 2 

arctan 

2 2pp l 



with 



-M W (E + iT w ) =p\ — ip 2 



(5) 



(6) 



For a complete calculation, one must add to eq. (§) the contributions from other final-state 
interconnection terms. 

For the energy region (|l]) one finds 



7T 

arctan x . 

2 



(7) 



In this region p\ and P2 have very transparent interpretatations: 

K T ~ — ~ — - 

pi r w 



(with p = p /M w ~ yjE/Myy) is the typical spatial separation between the diverging 
quasi-stable W's, while 

i i i 

(9) 



1 



1 



Pi Po \k\ 

is the characteristic distance of the Coulomb interaction between on-mass-shell W 
bosons. Thus the dominant contribution to 5c comes from values r^l/p , where 
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p = \J s/A — and r is the relative distance between the W's, as discussed in de- 
tail in ref. Hl4| . 5q depends crucially on the virtuality, and is controlled by the phase 
shift between the spatial oscillations corresponding to the actual momentum p and the 
characteristic momentum p 2 . Let us briefly discuss the main messages of an analysis of 
this momentum distribution. 

1. For stable W's the on-the-mass-shell Coulomb term would always induce a shift of 
the distribution towards lower momentum values. In this limit, the exact Coulomb 
result PB| is actually twice the first-order contribution (tt/2 in eq. (El)). 



For unstable W's 5q is always below the on-the-mass-shell value of tt/2 for p < \k\ 
and exceeds it for p > \k\. The transition occurs dominantly when 

\p~P2\ ~£>1 • (10) 

For region ([I]) this corresponds to \x\ £ 1. When p < p the effects of destructive 
interference accumulate in the course of integration over distances r up to WM 



p-po r w 

and, as a result, 5q gets screened. Correspondingly, when p > p , there is a con- 
structive interference in the same range of r values, leading to an increase of 5c- 
In the formal limit ^> \k\ 3> p one has 

C=po~a . 2P1 . (12) 
Contrary to the stable case, the Coulomb correction C does not blow up at (3 —>■ 



but reaches its maximal value of ay2Mw/Tw at E = 0. Below threshold, when 

\e\ r^, 

In region (|I|), relevant for LEP 2, Coulomb corrections are suppressed, 

The screening of the 1/p singularity for p p also takes place for higher-order 
Coulomb effects. This can be explicitly seen from the general formulae presented 



in ref. |[14|| . The higher-order Coulomb terms could in principle be more important 
numerically for differential distributions than for the total cross section; see however 
the discussion in the next section. 
4. 5c — > 7r, i.e. twice the leading-order on-the-mass-shell result, when M\y 3> p 3> 
However, when E ^> Tw and Mw ^ V 3> Po the Coulomb correction aix jp formally 
is much smaller than the on-the-mass-shell value aix / (2po). In region ([I]) the arctan 
modification of 5q in eq. (^) is an odd function of x, so after integration over p the 
total cross section agrees with the stable- W result. 
The instability effects induced by Coulomb interactions are expected to be reduced 
at relativistic energies because of the contributions of other final-state interaction mech- 
anisms. We will return to this issue below. 
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3 Numerical results 



In the following calculations, we assume a W mass M\y = 80.41 GeV, based on the recent 
CDF number |TJ5|], so as to simplify comparisons with the results of ref. [TO. The width 
is then T w = 2.1 GeV. Results are presented in the energy range 150-200 GeV, and for 
six discrete energies: 10 GeV and T\y below threshold, at threshold, Tw above threshold, 
at 175 GeV (typical LEP 2 running energy) and 30 GeV above threshold. Initial-state 
photon radiation is neglected, so as not to confuse the issue — in an idealized world 
this would be obtained by the process 77 — > W + W~. Other interfering graphs, such as 
e + e~ — > Z°Z° — > 4 fermions, are not considered. Only trivial loop corrections are included, 
such as the running of a and a s . The lowest-order cross section a® is therefore given by 



the Muta et al. expression pOR 



Several comparisons will be presented between alternative descriptions of Coulomb 
effects. This will help us assess the impact of QED interconnection, and to estimate 
uncertainties. 

• The simplest possible alternative is to have no Coulomb corrections at all, hence- 
forth designated "no Coulomb". Many studies in the past have been based on this 
approach, i.e. using a unmodified. 

• In the first-order "stable Coulomb" description, the two W bosons are considered 
as stable particles. Then a is modified by a simple factor 1 + a5c/P = 1 + na/2p. 

• The "unstable Coulomb" approach is a realistic standard for the effects of unstable 
W's. Again the multiplicative factor is 1 + aSc//3, but now 5q is given by eq. @. 

• Recently simple formulae were presented for the calculation of the Coulomb effects 
for unstable W's to second order ||21|| , "second-order Coulomb". The cross section 



is now given by <To multiplied by a factor |/(p, E)\ , with 

/(P, E) = 1 + ^ In D + #- f 1 ^ In l±fg_ , (15) 



4ip 16ipn Jo z 1 + z/D ' 

where D = (k + ip)/(n — ip) and k is defined in eq. (|6|). 

Fig. 1 compares the total e + e~ — > W + W~ cross section in the four scenarios above 
as a function of cm. energy. The differences are more readily visible in Fig. 2, where 
the cross sections have been normalized to the no-Coulomb one. The corrections reach 
a maximum at around threshold, of about 6%, and thereafter drop, asymptotically like 
1//3q. They also drop below threshold. This latter behaviour may be understood from 
the shape of the differential momentum distribution da /dp in the four scenarios above, 
Fig. 3. Note the very broad momentum spectra at energies below the threshold, where (at 
least) one of the W's is pushed significantly off the mass shell, and thus has a wide mass 
distribution. This depletion of the low-/? region thus reduces the Coulomb corrections 
below threshold. Above threshold the spectrum is peaked close to the on-mass-shell value 
p with a width decreasing like 5p ~ M w T w /p . 

As above, differences between the alternative Coulomb scenarios are more easily seen 
it the ratio of da /dp distributions with and without Coulomb effects included, Fig. 4. 
In particular note the characteristic destructive interference for p < po in the unstable 
Coulomb case. 

Even with this expanded scale, the difference between first- and second-order unstable 
Coulomb is very small. To further quantify effects, Fig. 5 gives the ratio of the second- 
order contribution to the first-order one. At energies below and around threshold, the 
effects are positive but only a few per cent, and vary slowly with p. Above threshold, the 
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difference between low and high p is noticeable on this scale. For p > p the effects are 
close to vanishing, < 1%, but become negative of 0(10%) for for p < p . In the limit 
p <C Po a simplified analytical expression is 



\f(p,E)f 



2aM w pi 



w 



2 Pl 



1 + 21n2 



pi 



vl 



(16) 



This gives a ratio indicated by a dashed line in Fig. 5. The sign of the second-order 
effect is such as to further reduce the next-to-vanishing first-order Coulomb corrections 
for p < p . We must therefore conclude that it will be impossible to observe any second- 
order Coulomb effects in the differential distributions at LEP 2. This alternative will 
largely be dropped from the continued discussion, although it is included in the figures. 

As noted earlier, the contributions of other intermediate-final and final-final inter- 
action mechanisms have been neglected in the above scenarios, either the calculation is 
first- or second-order. Driven by physical intuition we might try to estimate the possible 
impact of such screening corrections by using a simplified formulae containing an extra 
factor (1 — (3) 2 in front of the arctan term in eq. (||) of the first-order unstable Coulomb 
formulae, "dampened Coulomb". This is in agreement with the results of ref. [13], where 
an attempt is made to perform the quantitative studies of the instability effects in the 
relativistic region. However, it should be noted that the asymptotic behaviour does not 
constrain the form for LEP 2 energies, where (3 ~ 1/2, and thus the numerical difference 
between two asymptotically equivalent forms such as (1 — (3) 2 and (1 — /3 2 ) 2 (or (1 — (3) 2 
and (1 — A)) 2 ) is very large. The curves for this alternative at 175 and 190 GeV in Fig. 4 
should therefore be taken SIS db trial only. 

The ansatz above gives a Coulomb effect that is essentially constant as a function 
of p. Around po the sharp increase noted in the unstable Coulomb case is considerably 
dampened, since (1 — (3) 2 ~ (1 — 0.5) 2 = 0.25. The visual impression is of an even stronger 
dampening than just a factor 4, since the general 1//3 fall-off comes in addition and is 
common. Formally, the dampened Coulomb option can be seen as a mixture of (1 — (3) 2 
parts of the unstable Coulomb and 1 — (1 — f3) 2 of the stable Coulomb, which offers an 
alternative approach to understanding the role of the 1/(3. Since the fraction of the stable 
Coulomb vanishes like 2(3 — (3 2 for (3 — > 0, but the Coulomb factor increases like 1/(3, the 
flat behaviour for small p is to be expected. This should not be given too much emphasis, 
however, since we already noted that the whole approach is unreliable for the LEP 2 
energy region and in particular for momenta that deviate too much from po. However, in 
summary, it must be concluded that the dampened Coulomb scenario looks rather much 
like the no-Coulomb ones, except that the cross section is larger. So long as we are only 
concerned with distributions normalized per event, no separate figures need therefore be 
given for this alternative. This conclusion is not valid around or below threshold, but 
numerically is a sensible approximation from 170 GeV onwards. 

We now proceed to quantify how much event properties are changed by Coulomb 
effects. Fig. 6 shows the average W momentum as a function of cm. energy. At large 
energies this is increasing like (p) « \J s/A — M^, but around threshold one or both Ws 
are pushed off mass shell and so (p) stabilizes at a level around 20 GeV. The addition of 
Coulomb effects gives a shift of around 120 MeV in the threshold region, Fig. 7. The effect 
is here negative, i.e. the Coulomb factors enhance the low-momentum tail of the spectrum 
(cf. Fig. 4), giving a reduced (p). From 164 GeV onwards the destructive interference 
at low momenta leads to a net (p) shift upwards in the unstable Coulomb case. Lower 
momenta are always favoured for the (unrealistic) stable Coulomb alternative, so the 
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asymptotic convergence towards the no-Coulomb baseline is from opposite directions for 
a description in terms of stable or unstable W bosons. 

Fig. 8 shows the shift between the average W mass of events and the nominal M^f 
mass parameter (cf. ref. Around and below threshold this is a negative number, 

while it becomes positive at around 190 GeV. This is explained as a combination of phase 
space and cross section effects, as discussed in the introduction. Coulomb effects on the 
average W mass (m) are opposite in sign and about half as large as those on the average 
W momentum (p), Fig. 9. This is easily understood from the relation p 2 + m 2 ws/4 ==>• 
5m « —p5p/m, where 1/4 ^ (p)/m ^ 3/4 for the energy range considered here. The mass 
difference betweem first- and second-order Coulomb effects is seen to be below 3 MeV 
everywhere, i.e. completely negligible. 

Since each event contains two W's, it is possible to study separately the mass distri- 
bution of the lighter and the heavier W of each events. Figs. 10 and 11 show the Coulomb 
effects on the average light and heavy mass, min(m 1 ,m 2 ) and max (mi, m-z). Below and 
around threshold almost the whole W mass shift is found in the lighter W: since the 
lighter W here is off the mass shell, a further change of the mass means a smaller rel- 
ative change of the propagator value than a corresponding change for the heavier W. 
Mathematically, a W propagator p ~ l/((s — M^) 2 + M^T 2 ^) gives a relative change 
(l/p)|dp/ds| = 2\s-M^\/{{s-M^) 2 +M^T 2 N ) that is maximal around y/s ~ M w ±r w /2 
and falls off like l/\s — M^\ in either tail. At larger energies the mass shift is shared evenly 
between the two W's, which here both are about equally close to the nominal mass M\y. 

In the threshold region, there would then be two reasons to base a topological deter- 
mination of the W mass on the mass distribution of the heavier W of the event. Firstly, 
Coulomb corrections are small, and therefore also the associated uncertainty. Secondly, 
this mass distribution is significantly narrower that that of the lighter W or even that of 
the average mass (m heavy + m light )/2, Figs. 12, 13 and 14, so the statistical error would 
be smaller. However, there are also two reasons against such a study. Firstly, it would 
be a major experimental challenge to reconstruct two W masses per event, so the loss 
in accuracy from experimental effects would almost certainly outweight the above gain. 
Secondly, in this energy region a better accuracy can be obtained from a measurement of 
the total cross section ||. 

The root-mean-square width of the average W mass distribution is narrowest at around 
170 GeV, Fig. 14. Below this energy it increses significantly for the reasons discussed 
above; above 170 GeV there is a very slight increase from phase space effects. Coulomb 
effects narrow the mass distribution at around threshold, Fig. 15. Again this behaviour 
can be related to the width of the momentum distribution, Fig. 16, and to a narrowing of 
the width when Coulomb effects are included, Fig. 17. In general, however, the relation 
between a m — a m>n c and a p — cr p nC is not as straightforward as between (m) — (m nC ) 
and (p) — (pnc)- This is because the factor p in the 5m m —p5p/m relation suppresses 
the importance of the low-momentum tail on the mass distribution, while it enhances 
the importance of the high-momentum tail. Thus the stable-Coulomb scenario has a 
broader momentum distribution than the no-Coulomb one at most energies, but since 
this is due mainly to the enhancement of the low-momentum tail, the mass distribution 
may still be narrower. In the unstable-Coulomb approaches, the enhancement is in the 
high-momentum tail, and thus the momentum and mass width changes better follow suit. 

As we have noted earlier, experimental W mass methods are based on the p observable 
rather than on the mi and m 2 ones. That is, measured jet/lepton directions and energies 
are not used to reconstruct two masses per event. Instead the two W's are assumed 
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to have the same mass, and beam energy constraints are used to improve the detector 
energy resolution, so that most of the experimental information is related to angles and 
from there to the net momentum of jet pairs, i.e. to p. The translation from a set of p 
measurements to a W mass can be done in different ways, that are more or less sensitive 
to Coulomb effects. Compare Figs. 18 and 19; in the former the momentum p is first 



averaged over events and then translated to an average mass, (to) = Js/A— (p) 2 , in the 



latter the p's are first translated to to's and then averaged, (to) = (ys/4 — p 2 ). It turns 
out that the former recipe is preferrable in the threshold region, in that it gives a smaller 
spread between the alternatives, while there is no noticeable difference at higher energies. 

Other averaging schemes can be considered, where also events are giving relative 
weights to reflect their "usefulness" in constraining the W mass. The most direct ap- 
proach is to remove some fraction of the events in the wings of the momentum distribu- 
tions, where events are found only if one of the W's is significantly off the mass shell. For 
experimental rejection of background events, such an approach offers additional benefits. 
The importance of the wings for the momentum shift of Coulomb effects is shown in 
Fig. 20 for a few energies. What is plotted is 

dA(p> = d((p c ) - (P„c» ( J_ Ave _ J_ 4go\ 

dp dp \a c dp cr nC dp J 

where index "C" ("nC") denotes a scenario with (without) Coulomb effects. By construc- 
tion, the measure vanishes at (p) n cj this is visible at the lower energies but is somewhat 
hidden by the finite resolution at the highest energies. In other words, the shift of the 
average momentum is dominated by the behaviour in the wings at low energies. At 
higher energies, the momentum distribution is narrower and also the shift of the average 
momentum is dominated by a smaller region. This behaviour carries over to the mass 
shift (to) — (m n c). We study two alternatives, keeping only those events with momentum 
within 30 and 10 GeV, respectively, of the average. For simplicity, (p) has here been 
approximated by the expression (p) apP rox = \J{s/A — M^) 2 + 4M 2 v r 2 v , which is not par- 
ticularly accurate below threshold but good enough for this purpose. Comparing the mass 
shifts in Figs. 21 and 22 with those in Fig. 9, we note the reduced importance of Coulomb 
effects below and around threshold, while the reduction is rather modest at higher ener- 
gies. For instance, at a typical LEP 2 energy of 175 GeV, the difference between unstable 
Coulomb and no Coulomb decreases from —23 MeV without cut to —16 MeV with a cut 
\p ~~ (p)approx| < 10 GeV. The statistics will drop if the momentum window is reduced 
further, so it appears difficult to do much better. 

We have compared several alternative Coulomb effect scenarios in this paper. Thereby 
we hope to learn about the possible effects also of other (not yet calculated) QED effects, 
like intermediate-final interference. Whereas a scenario such as the no-Coulomb one 
obviously is excluded, we have above indicated how something akin to it might arise as 
an approximate description of the normalized momentum distribution. Such ambiguities 
would be resolved if the experimental data themselves could be used to extract the form 
of QED corrections. One specific idea would be to divide the experimental momentum 
spectrum by the no-Coulomb rate, to obtain a distribution akin to what is shown in Fig. 4. 
If the experimental distribution would agree with contemporaneous calculations, theory 
could be used to extract the W mass. If not, the experimentally defined correction factors 
could still be used for this purpose. 

There are two major problems with such a thinking, however. Firstly, so long as we 
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do not know the W mass we also do not know which no-Coulomb curve to compare with. 
The sharp increase in the unstable Coulomb factor at around po for the higher energies 
in Fig. 4 is there since we compare scenarios with the same M^. If instead (p) is made to 
agree, as would be the realistic experimental procedure, differences between scenarios are 
not as spectacular. Still, even with the same (p), the shape of the momentum spectrum 
would be different, so with hard work it should be possible to find which descriptions 
work and which do not. However, here enters the second problem, namely that of limited 
statistics. A LEP 2 experiment will at most have of order 10,000 events to base the 
analysis on. It is difficult to see how this would be enough to observe any differences, 
given that we are speaking of Coulomb correction factors of a few per cent. 

This prejudice is borne out by detailed studies. For instance, in Fig. 23 we have com- 
pared 100 "LEP 2 experiments" at 175 GeV, each of 10000 events. In the stable-Coulomb 
scenario the Mw parameter has been shifted down by 13 MeV relative to the no-Coulomb 
one, and in the unstable-Coulomb one up by 22 MeV (cf. Fig. 9), so that (p) should be 
the same. In fact, the limited statistics leads to a non-negligible spread of (p) values. Also 
higher moments, in Fig. 23 \j ({p — (p)) n ), n = 2,3, 4, show a significant statistical spread, 
with no obvious separation between the no-, stable- and unstable-Coulomb scenarios. For 
a given "LEP 2 experiments" it is impossible to tell which distribution it is drawn from. 
Correlations between the moments again do not distinguish (not shown). 

One therefore must conclude that — unless something truly spectacular happens - 
experimental input can not be used to pin down QED interconnection effects. It is 
necessary to have the theory under control to the required accuracy. 



4 Conclusions 

One of the most critical single observables for LEP 2 physics is the W mass. Had there 
been no interconnection effects between the W + and the W~, in principle their respective 
four-momenta could be reconstructed and squared to give the W masses. The shifts in 



the observable W mass distributions could come from a number of sources || g, [lOj that we 
have not attempted to address here. Our main concern has been to estimate the W mass 
shift caused by distortions of the W momentum spectrum from instability effects, as em- 
bodied in the QED radiative corrections to the differential e + e~ — > W + W~ cross section. 
In the threshold region we can exploit the formulae for Coulomb corrections 0, [I4]| , which 
here are responsible for the dominant off-shellness-induced phenomena. For complete pre- 
dictions in the whole LEP 2 energy range, one needs to incorporate the contributions of 
other QED interconnection mechanisms, e.g. those resulting from intermediate-final and 
final-final photon interference || [13] . As far as we are aware, the complete analytical cal- 
culation has never been performed for the whole set of electroweak corrections to unstable 
W production and decay. (An attempt to estimate QED interconnection phenomena in 
the relativistic region can be found in ref. [|i~3||.) 

While the complete calculation is missing, one may attempt to compare some different 
scenarios that give an impression of the size of effects. Our current "best" description 
is the all-orders Coulomb description in ref. |21fl . One conclusion of this paper is that 
second-order (and, by implication, all higher-order) Coulomb effects are practically neg- 
ligible. The first-order ("unstable") Coulomb description therefore is a realistic baseline. 
Other intermediate-final and final-final interference terms are expected to modify this be- 
haviour. Specifically, one could expect the rapid variation of the Coulomb factor around 
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Po (the on-mass-shell W momentum) to be substantially dampened. In fact, a brute-force 



application of a proposed asymptotic dampening factor [13| to the LEP 2 energy region 
would imply that a momentum-dependent Coulomb factor can be replaced by an almost 
momentum-independent average. Therefore the no-Coulomb scenario could be a realistic 
alternative for reconnection effects on kinematical distributions. The final scenario, with 
Coulomb effects evaluated for stable W's, is less realistic and has been included mainly 
for completeness. 

At LEP 2 energies, the no- and unstable-Coulomb scenarios differ by about 20 MeV 
on the observable (m), given a common theory input Mw- This is somewhat below the 
statistical accuracy, but comparable with other potential sources of systematic error, and 
so not negligible. It is therefore interesting to see whether the uncertainty can be reduced 
by a clever choice of experimental procedure. The conclusions of this paper are not too 
optimistic on this count. The limited statistics will preclude the testing of the theory 
scenarios from data itself. One will therefore be totally dependent on a full-scale theory 
calculation to reach any definite conclusions. However, different statistical treatments of 
the data can increase or decrease the 20 MeV number above by some amount: in the 
paper we have given the examples of cuts on the momentum range analyzed and choice 
of statistical averaging procedure. 

The QED final-state interaction could induce some systematic effects in other W- 
mass measurements as well, for instance in pp collider experiments. Of particular interest 
here is the subprocess qg — > Wq' with W — > Ivg. Collider experiments normally rely on 
the equivalent process for Z° production, qg — ► Zq with Z — > £ + £~ , to calibrate the W 
mass scale. Non-universal interference effects are not included in such a procedure, e.g. 
a charged (W*) versus a chargeless (Z°) intermediate state. The effects of real-photon 



emission have been compared ]2"2"|| , but virtual corrections remain to be studied. It could 
well turn out that a corresponding theory uncertainty of order 20 MeV exists in this 
process. Within current experimental errors this would be negligible, but it could become 
relevant for future high-precision measurements. 
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;ure Captions 

1 Total e + e~ — > W + W~cross section a as a function of the cm. energy. Dashed- 
dotted curve is no Coulomb, dashed stable Coulomb, full unstable Coulomb 
and dotted second-order Coulomb. 

2 Total cross section relative to the no-Coulomb one, a /(Tq. Notation as in Fig. 1. 

3 The differential cross section da /dp in pb/GeV. Notation as in Fig. 1. 

4 Ratio of the differential cross sections with and without Coulomb effects in- 
cluded, (da"/dp)c/(dcr/dp) n c. Dashed curve is stable Coulomb, full unstable 
Coulomb and dotted second-order Coulomb. At the two highest energies, 
dashed-dotted curve is dampened Coulomb. 

5 Ratio of second- to first-order Coulomb corrections, (\f(p,E)\ 2 — 1 — 
oiSc/P)/(oi5c/P), as a function of p. The dashed line at small p indicates 
the p — > limit given in eq. ( |T6| ) . 

6 Average W momentum, (p). Notation as in Fig. 1. 

7 Shift of the average W momentum relative to the no-Coulomb value, 
(p) ~ (Pnc)- Notation as in Fig. 1. 

8 Shift of the average W mass relative to the nominal one, (m) — M^. Notation 
as in Fig. 1. 

9 Shift of the average W mass relative to the no-Coulomb value, (m) — (m nC ). 
Notation as in Fig. 1. 

10 Shift of the lighter W mass relative to the no-Coulomb value, ((m) — (m nC )) light . 
Notation as in Fig. 1. 

11 Shift of the heavier W mass relative to the no-Coulomb value, 
((m) — (m nC )) heavy . Notation as in Fig. 1. 

12 Root-mean-square width of the lighter W mass distribution, <7 mi ii g ht- Notation 
as in Fig. 1. 

13 Root-mean-square width of the heavier W mass distribution, cr mi h ea vy Nota- 
tion as in Fig. 1. 

14 Root-mean-square width of the average W mass distribution, a m . Notation as 
in Fig. 1. 

15 Shift of the width of the average W mass distribution relative to the no- 
Coulomb value a m — a m ,nc- Notation as in Fig. 1. 

16 Root-mean-square width of the W momentum distribution, a p . Notation as in 
Fig. 1. 

17 Shift of the width of the W momentum distribution relative to the no-Coulomb 
value, <7p — cr Pi nc. Notation as in Fig. 1. 
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Fig. 18 Shift of the average W mass relative to the no-Coulomb value, (m) — (m n c). 

Here the average momentum (p) has been converted into an average mass 
(m) = J s/4 — (p) 2 . Notation as in Fig. 1. 

Fig. 19 Shift of the average W mass relative to the no-Coulomb value, (m) — (m n c)- 
Here the momentum p in each event has been converted into an effective mass 



m = ys/4 — p 2 , and this mass has been averaged over events. Notation as in 
Fig. 1. 



Fig. 20 The buildup of the total shift in (p) as a function of p, as given in eq. ([TJ 
The result is given in relation to the no-Coulomb alternative. Dashed curve is 
stable Coulomb, full unstable Coulomb and dotted second-order Coulomb. 

Fig. 21 Shift of the average W mass relative to the no-Coulomb value, (m) — (m n c)- 
Here only events with |p — (p) ap prox| < 30 GeV have been included. Notation 
as in Fig. 1. 

Fig. 22 Shift of the average W mass relative to the no-Coulomb value, (m) — (m n c)- 
Here only events with |p — (p) ap prox| < 10 GeV have been included. Notation 
as in Fig. 1. 

Fig. 23 The distribution of some moments of the momentum distribution in 100 
"LEP 2 experiments" of 10000 events each at 175 GeV. Dashed-dotted his- 
togram is no Coulomb, dashed stable Coulomb, and full unstable Coulomb. 
The M\v mass parameter has been adjusted to give the same (p) in the three 
scenarios. 
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